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1.  Introduction 

In  numerical  weather  prediction  problems  it  is  customary  to 
apply  boundary  conditions  of  the  form  w  =  0,  w  =  0,  or  5  ®  0  or 
a  modification  of  either  of  these  at  the  upper  level.  Examples 
are  Kasahara  and  Washington  (1971),  Somerville  et  al.  (1974), 
Ramanathan  and  Grose  (1977),  Staniforth  and  Daley  (1977)  and 
Donahue  (1980).  Hie  purpose  of  this  report  is  to  highlight  some 
difficulties  of  these  boundary  conditions  and  to  suggest  a 
possible  alternative. 

It  is  convenient  to  consider  the  boundary  conditions  that 
have  been  applied  under  two  separate  headings: 

(a)  an  infinite  atmosphere  where  boundary  conditions  are 
applied  at  2  =  •  (identical  to  P  »  0),  and 

(b)  a  finite  atmosphere. 

(a)  Top  boundary  condition  for  an  infinite  atmosphere 

For  the  discussion  to  follow  it  is  worthwhile  to  note  the 
difference  between  a  true  boundary  condition  and  a  statement  of 
fact.  As  an  illustration  we  consider  the  solution  of  the 
equation 

^2+  28  |2  +  .Q  =  0  (X). 

d  Z 

If  both  a  and  8  are  real  and  positive,  a  statement  like  "Q 
tends  to  zero  as  Z  tends  to  infinity"  is  not  a  boundary 
condition;  rather  it  is  a  statement  of  fact.-  This  is  because 
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there  is  no  way  of  combining  the  two  different  possible  solutions 
to  Bq.(l)  without  necessarily  making  Q  tend  to  zero  as  Z  tends  to 
infinity. 

On  the  other  hand,  if  a  were  negative,  the  statement  would 
be  a  boundary  condition.  For  in  this  latter  case  one  of  the  two 
possible  solutions  increases  very  rapidly  as  Z  tends  to  infinity 
while  the  other  tends  to  zero.  The  statement  "Q  tends  to  zero  as 
Z  tends  to  infinity"  therefore  discriminates  between  the  two 
solutions. 

Now  let  us  consider  the  linearized  equations  of  motion  in 
the  P-system  of  co-ordinates; 


3  u 

7E  " 


fv 


0 


3  v 

TE 


+  fu  + 


dip1 

W~ 


o 


3  T' 
3 1 


+ 


d0  T« 
o  o 

3P“  W~ 
o 


=  0 


=  0 


(2). 


3u  ^  3v  ,  3w 
7x  +  3y  +  TP 


=  0 


Here,  u,  v  and  u  are  the  two  horizontal  components  of  velocity 
and  DP/Dt  respectively,  P  is  pressure,  41  and  T  are  geopotential 
and  temperature,  f  is  Coriolis  parameter  and  R  the  gas  constant. 

The  primed  quantities  are  perturbations  over  those  with 
subscript  zero  such  that  for  example  T  «  TQ  +  T' ,  | T* |  «  TQ. 

Some  algebraic  manipulations  yield 
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P2  ^  ) 

H2N2  3P 


*  liLZi  i!n 

N2H2  3  P 


where  a  solution  of  form 
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Tt 


iv 


is  assumed  and 


*’  = 


w 


z 

n 


“n* 


(3) 


Each  of  the  i|>n  and  u>n  are  the  spectral  Hough  components 
of  and  o)  respectively  and  hn  is  the  corresponding  equivalent 
depth  defined  in  Iaplace's  tidal  theory  (see  for  example  Wilkes, 
1949);  H  is  pressure  scale  height  and  N,  the  Brunt-Vaisala 
frequency,  is  given  by 


N2H2 

P2 


RT 

ae_ 

O 

o 

peo 

dP 

(4). 


Consider  an  isothermal  atmosphere.  The  solutions  to  Eqs. 
(3)  give 


m 


n 


m  =  1/2 


i  (1/4  -1^- 


Y-l  H 

>v 


(5) 


where  y  is  the  ratio  of  the  two  principal  specific  heats  of  air. 

Now  mesoscale  phenomena  and  a  large  number  of  synoptic  scale 
motions  have  positive  equivalent  depth.  From  Eq.  (5)  it  is  clear 
that  if  hn  is  positive,  the  real  parts  of  the  two  possible  values 
of  m  are  necessarily  positive.  That  is  to  say  that  as  P  tends  to 


zero  (■>  tends  to  zero  for  all  possible  solutions. 

We  could  have  arrived  at  this  conclusion  intuitively  as 

follows.  Following  the  argument  of  Staniforth  and  Daley  (1977) 

that  T  remains  bounded  as  P  tends  to  zero,  and  noting  from  Eq. 

T  de 

(4)  that  g—  -r=-  tends  to  infinity  as  P  tends  to  zero,  we  can  see 
°o 

immediately  from  the  third  equation  in  Eq.  (2)  that  u  is 
necessarily  zero  at  P  =  0. 

Thus  we  see  that  the  statement  "u>  =  0  at  P  =  0"  is  in  fact 
not  a  boundary  condition  for  many  of  the  possible  motions  in  the 
atmosphere.  Rather  for  such  motions  the  statement  is  only  a 
property  of  the  equations  of  motion. 

Similarly  if  a  =  P/Pg  where  Pg  is  surface  pressure,  the 

statement  "a  =  0  at  a  =  0"  is  not  a  boundary  condition  for 

mesoscale  and  other  motions  of  positive  equivalent  depths.  This 
can  be  shown  as  follows.  Note  that 

o  =  (d/P  -  P  f»  /p£  . 
s  s  s 

Since  in  =  0  at  P  =  0,  i.e.  at  a  =  0  ,  the  first  part  of  R.H.S.  of 

this  equation  is  zero  at  o  =  0.  The  second  part  is  also  zero 

because  P  =  0  at  a  =  0.  Therefore  a  is  necessarily  zero 
at  o  =  0;  and  so  once  again  the  statement  o  *  0  at  o  *  0  cannot 
be  a  boundary  condition. 

(b)  Top  boundary  condition  at  P  >  0 

Let  us  solve  Eq.  (3)  completely.  As  an  illustration,  apply 


the  boundary  conditions 
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u.  at  P  = 
s 


s' 


and  au  +  =  0  at  P 


where  Pg  and  PT  are  respectively  the  surface  and  top  pressures  of 
the  model  atmosphere  and  a  and  b  are  some  constants  to  make  the 
top  boundary  condition  general.  The  solution  of  Eq.  (3)  becomes 


(i> 


n 


p  M  A<p/p  ix  -  bip/pt)-x 

‘V'P-™  l - X - =T' 

s  *s  a<ps/pt)a  -  b(ps/pt) 


(6) 


where  A  =  aPT+  -X)b,  B  =  aPT  +  +  X)b,  X2  *  V4  -  Ili  JL. 

T  n 


We  can  discuss  Eq.  (6)  under  the  two  special  cases  of  X  real 
and  X  imaginary. 


Case  (a):  X  real 

This  case  is  relevant  for  large  equivalent  depths.  If  the 
top  is  high,  i.e.,  PT  very  small,  Eq.  (6)  at  lower  levels  become 

•„  *  (p-fc  (|-)X  (7) 

s  s 

which  is  independent  of  the  top  boundary  condition. 

Thus  provided  the  level  (at  which  the  top  boundary  condition 
is  applied)  is  high  enough  and  Eq.  (3)  is  not  oscillatory  in  the 
vertical,  the  particular  boundary  condition  and  the  level  at 
which  it  is  applied  have  little  or  no  effect  on  the  solution  to 
the  equation. 


This  is  the  case  for  equivalent  depths  such 
Y  —  1  X 

that  hn  <  4  -ly-  H.  in  this  case  (P/Pq,)  is  not  only  a  reciprocal 
of  (P/PT)“*  but  also,  is  a  complex  conjugate  of  the  other.  So 
reducing  Pq.  does  not  make  one  solution  more  important  than  the 
other  -  rather  both  solutions  go  through  series  of  maxima  and 
minima  as  PT  is  altered.  The  overall  effect  is  that  o>n  at  a 
particular  P  is  critically  dependent  on  Pq,  and  on  a  and  b.  That 
is,  the  solution  depends  very  strongly  on  the  boundary  condition 
and  the  level  at  which  it  is  applied.  This  is  clearly 
undesirable  as  it  makes  the  final  solution  very  sensitive  to  the 
parameter  (boundary  condition)  of  which  we  are  uncertain. 

It  may  also  be  noted  from  Eq.  (3)  that  the  boundary 
condition  u  =  0  implies  that  d$'/3P  =  0.  Apart  from  W  =  0  which 
may  be  interpreted  as  a  rigid  lid  representing  the  inversion  at 
the  tropopause,  it  is  difficult  to  think  of  some  physical 
processes  represented  by  other  boundary  conditions  like  u  =  0  or 
a^'/SP  =  0  or  o  =  0. 

ihis  difficulty  in  choosing  the  upper  level  boundary 
condition  is  common  in  geophysical  problems.  In  the  case  of  the 
atmospheric  diurnal  tides,  the  solution  was  to  choose  a  condition 
that  was  consistent  with  the  effect  of  damping  due  to  diffusion 
of  heat  and  momentum;  see  for  example  Pekeris  (1937),  Giwa 
(1967).  This  procedure  is  reasonable  since  whenever  there  is 
some  motion  there  is  necessarily  some  diffusion.  Therefore  in 
the  next  section  we  shall  introduce  vertical  eddy  diffusion  of 
the  horizontal  components  of  velocity  and  heat  into  the 
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linearized  equations  of  motion. 


2.  Inclusion  of  eddy  diffusion  in  the  linearized  primitive 
equations  of  motion 


(a )  Governing  equations 

In  order  to  simplify  the  derivation  to  follow  we  shall  in 
the  first  instance  consider  meso-scale  motions;  that  is  we  shall 
neglect  Coriolis  acceleration  in  comparison  to  3/3t  terms.  It 
will  be  shown  later  that  the  conclusions  arrived  at  will  be 
equally  applicable  if  the  Coriolis  terms  are  included.  Indeed  it 
can  be  deduced  that  the  effect  of  the  Coriolis  terms  is  to  modify 
the  horizontal  constant  of  separation  which  is  often  represented 
by  the  equivalent  depth. 

Using  standard  notation  the  linearized  primitive  equations 
of  motion  in  the  "P"  system  are 


3  u  .  3  + ' 
"ST  +  Tx~ 


=  (npQ  fp) 


3v  .  3*' 
TT  +  Ty“ 

3*’  .  RT' 
TP~  +  T 


=  g2  TP  (npo  Ip) 


3  T'  d0o  T0 

Tt~  +  “  3p"  e“ 
o 


)2  3  fv  3T'i 

b  I?  (Kl>o  Tp-] 


3 U  ,  3 V  ,  3u 
Tx  +  7y  +  TP 
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where  g  is  acceleration  due  to  gravity,  Cp  is  specific  heat  at 
constant  pressure,  n  and  K  are  the  eddy  coefficients  of  viscosity 
and  conductivity  and  the  remaining  symbols  are  as  defined  in  Eq. 
(2). 

From  Eqs.  (8)  and  (9)  we  have 


( np. 


3 

3P 


>  w3u  3V, 
J  ;t3x  3yf 


+  V2*'  =  0 


V2 


3  2  +  3  2 

3x2  3y2 


(13). 


Differentiating  (13)  with  respect  to  P  and  using  (12)  to 
eliminate  we  have 


3 

TV 


[{ 


3 

TV 


~9Z 


3 

TP 


(  npo  Ip) 1 


3uj  -{ 

W 


V 2  T’  =  0 


(14) 


Combining  Eqs.  (11)  and  (14)  to  eliminate  T'  we  have 


r  3  g2  3  r  v 

'Ft  ”  C  Tp(Kp 
P 


3  1  IfZ  3 
o  Tpj '  lr  IP 


3  2  w 

TETp 


-  „2i 


tnp„  -  (15) 


TP'-MKo 


3P" 


T  de 

o  o 

e  3p~ 

o 


V2u>  =  0 


Now  the  purpose  of  this  discussion  is  to  investigate  the 
modification  of  the  basic  equations  by  diffusion.  Therefore  we 
shall  assume  that  the  diffusion  terms  are  an  order  of  magnitude 
less  than  the  other  terms.  In  particular  we  shall  neglect 
products  of  the  diffusion  terms.  Also  we  shall  consider  an 
isothermal  atmosphere  in  order  to  obtain  analytic  solutions  in 
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terms  of  simple  functions.  It  may  be  noted  that  upper  boundary 
conditions  are  applied  in  the  region  of  the  tropopause  which  is 
very  nearly  isothermal. 

With  these  approximations  and  assuming  a  temporal  variation 
of  the  form  exp(ivt)  and  noting  that  in  an  isothermal  atmosphere 


de  T  _  y-1  gH 

cTP  TT  y  RP  ' 


Eq.  (15)  becomes 


+  |p(KP0  U&L))  ♦  1^0  0)1] 


-  V2u>  =  ° 


If  for  the  sake  of  obtaining  an  analytic  solution  we  make  the 
reasonable  assumption  that  the  kinematic  coefficients  of 
diffusion  are  constant  with  height,  i.e.,  K/pQ  =  K*/p*  =  constant 
and  n/pQ  =  n*/p*  •  constant,  and  noting  that  the  equivalent  depth 
h  is  given  by 


g?2  =  -v2/hn  , 


we  have 


'p  H2 


p  3  ( p2  3  pa 

~  tpp- 


pz  a' 


**  h2  ap2  ap2 


*  a)  * 


Y-1  H 

Y  TT  n 
n 


which  can  be  simplified  to  give 

3  2u  3  **  w  3 3  to  32w  .  __ 

p2 - £  +  iN  (p* - -  +  4P3 - -  +  2P2 - 2)  +  Hi  *L  u  m  0 

*  »2  an1*  ao3  Y  H  n 
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where  the  non  dimensional  constant  ND,  representing  the  effect  of 
diffusion,  is  given  by 


1  rK*  n*x 

N°  ’  ^  • 


(b)  Solution  of  the  equation  and  application  to  upper  level 
boundary  condition 

Eq.  (17)  is  homogeneous  in  P  and  hence  we  can  seek  a 
solution  of  form 

u>  «  PM  (18) 

n 

When  this  is  substituted  into  Eq.  (17)  we  have 

M(M-l)  +  iNDM2  (M— 1 ) 2  +  £  =  0  (19). 

The  discussion  to  follow  is  perhaps  best  done  in  terms  of 
the  Z  system  of  coordinates.  ALso  it  is  convenient  to  define  a 
new  quantity 

Q  *  “/Pq  ( 20 ) 

Since  P  *  P  e”Z/^H  (P_-mean  surface  pressure),  substituting  Eq. 

OS  ® 

(18)  into  Eq.  (20)  gives 

0  «  eXZ/H;  X  -  1-M  (21) 


so  Eg.  (19)  in  terms  of  X  gives 


X(X-l)  +  iNDX2(X-l)2  +  £  =  0  (22). 

Eq.  (22)  is  a  fourth  order  equation  and  hence  has  four 
possible  independent  solutions,  two  of  which  do  not  exist  if  the 
damping  coefficients  are  set  to  zero.  Henceforth  we  shall  refer 
to  these  two  solutions  as  the  diffusion  terms.  The  other  two 
solutions  will  exist  even  when  the  coefficients  of  diffusion  are 
zero  and  they  are  the  solutions  which  have  always  represented  the 
normal  modes  of  atmospheric  motion.  The  purpose  of  this 
discussion  is  to  see  how  these  other  two  solutions  are  modified 
as  a  result  of  the  inclusion  of  diffusion. 

Before  proceeding  with  the  discussion  it  is  worthwhile  to 
note  some  salient  points. 

(1)  The  total  kinetic  energy  in  a  column  of  air  from  the 
bottom  to  the  top  of  the  atmosphere  (at  Z  =  ZT)  is  proportional 
to 

^T  n 

/  (tp  dZ  where  Q  is  defined  in  Eq.  (20).  Since  the 
o 

atmosphere  is  of  finite  mass,  the  total  energy  must  be  finite 

even  as  ZT  tends  to  infinity.  This  means  that  Q  should  increase 

1a 

less  rapidly  with  Z  than  decreases.  From  the  expression  for  0 
in  Eq.  (21)  and  since  pQ  «  e“Z^H,  the  finiteness  of  kinetic 
energy  means  that  the  real  part  of  *  should  be  less  than  • 

(2)  A  wave  that  is  generated  as  a  result  of  reflection  from 
the  top  should  have  its  amplitude  decreasing  more  rapidly  (as  it 
propagates  downwards  through  the  viscous  atmosphere)  than  if 


there  were  no  diffusion.  The  obvious  implication  of  this  is  that 
the  amplitude  of  a  wave  reflected  from  above  should  increase 
faster  in  the  positive  vertical  direction  than  the  fundamental 
wave  whose  source  is  the  energy  from  the  lower  levels. 

Table  1  gives  solutions  of  Eq.  (22)  for  different  values  of 
H/h  and  ND.  The  first  value  for  each  X  is  the  real  part  and  the 
second  represents  the  imaginary  part.  It  should  be  noted  that 
since  Eq.  (22)  is  fourth  order  in  X  there  are  four  different 
X's  for  each  solution  except  the  special  case  ND  =  0,  which  has 
two  different  X*s.  From  the  table  it  can  easily  be  seen  that 
the  first  set  of  values  for  X  in  each  box  are  physically 
impossible  since  these  will  result  in  unbounded  energy  as  Zrp 
increases.  This  set  is  in  fact  one  including  the  diffusion 
terms.  We  can  easily  identify  the  solutions  involving  diffusion 
terms  from  an  examination  of  Eq.  (22).  As  ND  is  made  smaller  and 
smaller  the  corresponding  X  for  the  diffusion  terms  increases  in 
magnitude  very  rapidly  -  whereas  X  for  the  other  two  terms  tends 
to  a  solution  of  the  equation 


‘2  -  »  +  ^  K 


0. 


Thus  we  can  see  that  the  second  set  of  values  in  Table  1  also 
belong  to  the  diffusion  term. 

Let  us  for  convenience  call  the  remaining  set  of  values 
Xft  and  X gj  XA  being  the  3rd  set  and  XB  the  fourth  set.  A 
statement  of  our  problem  is  to  see  which  of  XA  and  Xfi  is 
physically  realizable  in  an  atmosphere  with  an  open  top.  From 
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Table  1  we  see  that  the  real  part  of  *A  increases  as  the 
coefficient  of  diffusion  increases.  Thus  this  wave  has  its 
source  of  energy  from  above  (i.e.  reflected  wave)  and  in 
propagating  downwards  it  damps  out  more  readily  as  the 
coefficients  of  diffusion  are  increased. 

For  Ag  the  real  part  decreases  as  the  diffusion  coefficients 
increase.  Hence  this  is  a  wave  coming  from  the  lower  levels. 

This  is  in  fact  the  only  wave  that  satisfies  the  following 
conditions: 

(1)  It  does  not  vanish  when  the  coefficients  of  diffusion 
are  zero; 

(2)  its  existence  is  not  a  consequence  of  reflection  from  the 
top. 

Ihus  the  suggested  boundary  condition  at  the  top  is  that 
/pQ  Q  shall  behave  as  exp  i(vt  +  kz).  This  means  that  at  the 
top  the  phase  velocity  should  be  directed  downwards.  It  may 
appear  paradoxical  that  the  phase  velocity  is  directed  downwards 
when  the  source  of  energy  is  from  below.  As  shown  by  Wilkes 
(1949)  such  a  system  is  in  fact  consistent  with  energy 
propagating  upwards. 

We  note  that  since  the  solution  does  not  contain  either  PT 
or  Z.p,  the  boundary  condition  is  independent  of  the  level  at 
which  it  is  applied. 
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Table  1:  Complex  solutions  for  X  (Eq.  22)  for  various  values  of  and  H/h. 


nd 

0 

0. 

01 

0. 

1 

1. 

0 

10 

•  0 

H/h 

7.58 

7.06 

2.76 

2.21 

1.30 

0.62 

1.01 

0.10 

-6.58 

-7.06 

-1.76 

-2.21 

-0.30 

-0.62 

-0.01 

-0.10 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

T.  58 

■tKEfli 

■EKU 

171 1 

1.31 

0.62 

MTTT  ■ 

■iWHl 

-6.58 

-7.06 

-1.77 

-2.21 

-0.31 

-0.62 

-0.03 

-0.10 

0.10 

0.97 

0.00 

0.97 

0.00 

0.97 

0.00 

0.97 

0.00 

0.97 

-0.01 

0.03 

0.00 

0.03 

0.00 

0.03 

0.00 

0.03 

0.00 

0.03 

0.01 

7.58 

7.06 

2.77 

2.20 

1.33 

0.61 

1.06 

o.ii 

-6.58 

-7.06 

-1.77 

-2.20 

-0.33 

-0.61 

-0.06 

-0.11 

0.25 

0.92 

0.00 

0.92 

0.00 

0.92 

-0.00 

0.92 

-0.01 

0.95 

-0.02 

0.08 

0.00 

0.08 

0.00 

0.08 

0.00 

0.08 

0.01 

0.05 

0.02 

7.58 

7.06 

2.78 

2.19 

1.36 

0.59 

1.09 

0.12 

-6.58 

-7.06 

-1.78 

-2.19 

-0.36 

-0.59 

-0.09 

-0.12 

0.50 

0.83 

0.00 

0.83 

0.00 

0.83 

0.00 

0.84 

-0.03 

0.92 

-0.05 

0.17 

0.00 

0.17 

0.00 

0.17 

0.00 

0.16 

0.03 

0.08 

0.05 

7.59 

7.06 

2.80 

2.18 

1.42 

0.58 

1.12 

0 . 14 

-6.59 

-7.05 

-1.80 

-2.18 

-0.42 

-0.58 

-0.12 

-0.14 

1.00 

0.50 

-0.19 

0.50 

-0.19 

0.52 

-0.19 

0.67 

-0.18 

0.88 

-0.10 

0.50 

0.19 

0.50 

0.19 

0.48 

0.19 

0.33 

0.18 

0.12 

0.10 

7.61 

7.04 

2.85 

2.14 

1.54 

0.58 

1.18 

0.18 

-6.61 

-7.04 

-1.85 

-2.14 

-0.54 

-0.58 

-0.18 

-0.18 

2.50 

0.50 

-0.68 

0.50 

-0.68 

0.54 

-0.68 

0.69 

-0.54 

0.83 

-0.21 

0.50 

0.68 

0.50 

0.68 

0.46 

0.68 

0.31 

0.54 

0.17 

0.21 

7.63 

7.01 

2.94 

2.08 

1.68 

0.60 

1.25 

0.2l 

-6.63 

-7.01 

-1.94 

-2.08 

-0.68 

-0.60 

-0.25 

-0.21 

5.00 

0.50 

-0.09 

0.51 

-1.09 

0.59 

-1.07 

0.77 

-0.77 

0.82 

-0.34 

0.50 

1.09 

0.49 

1.09 

0.41 

1.07 

0.23 

0.77 

0.18 

0.34 

7.68 

6.97 

3.12 

2.02 

1785" 

6.65 

1.32 

0.26 

-6.68 

-6.97 

-2.12 

-2.02 

-0.85 

-0.65 

-0.33 

-0.26 

10.00 

0.50 

-1.62 

0.52 

-1.61 

0.70 

-1.53 

0.87 

-1.02 

0.84 

-0.49 

0.50 

1.62 

0.47 

1.61 

0.30 

1.53 

0.13 

1.02 

0.16 

0.49 

S77T 

6 . 46 

"  4.54 

2 .  l6 

'  2 . 72 

or.  97 

— rrrr 

O'.  46  j 

-7.71 

-6.46 

-3.54 

-2.16 

-1.72 

-0.97 

-0.79 

-0.48 

100.00 

0.50 

-5.32 

1.10 

-5.06 

1.57 

-3.51 

1.31 

-2.04 

1.02 

-1.10 

0.50 

5.32 

-0.10 

5.06 

-0.57 

3.51 

-0.31 

2.04 

-0.02 

1.10  7 
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3.  Comparison  with  similar  work  and  conclusion 

Lilly  and  Kennedy  (1973)  observed  that  waves  with  small  wave 
numbers  tend  to  propagate  through  the  troposphere  with  little 
reflection  of  energy.  This  led  Klemp  and  Lilly  (1975)  to  apply  a 
boundary  condition,  which  allowed  energy  to  propagate  upwards. 

Thus  the  observation  of  Lilly  and  Kennedy  agreed  with  our 
theoretical  discussion  here  and  an  old  suggestion  of  Weekes  and 
Wilkes  (1947)  in  which  they  pointed  out  that  since  most  of  the 
tidal  energy  is  supplied  at  lower  levels,  the  direction  of  energy 
flow  at  ZT  should  be  positive. 

Lindzen  (1974)  applied  a  radiation  condition  at  the  top  of 
the  atmosphere.  He  applied  a  condition  which  made  wave  velocity 
positive  in  the  direction  of  increasing  Z.  The  basis  for 
preference  of  positive  wave  velocity  as  opposed  to  negative  wave 
velocity  was  not  discussed.  On  the  other  hand,  as  discussed 
above,  Wilkes  (1949)  showed  that  a  negative  wave  velocity  was 
consistent  with  energy  propagating  upwards. 

Using  the  Sommerfeld  radiation  condition,  Orlanski  (1975) 
proposed  a  boundary  condition  by  assuming  that  each  of  the 
variables  satisfies  a  wave  equation  of  form 


1  ii  = 
cat  TZ 


at  the  top  boundary.  Hie  value  of  C  was  calculated  from  previous 
values  of  t  at  the  boundary  and  one  spatial  grid  away  from  the 
boundary.  Orlanski  showed  that  by  this  procedure,  if  the  problem 
was  linear  and  consisted  of  a  single  wave,  there  was  no  reflected 


wave.  If  however  there  is  more  than  one  wave,  as  is  always  the 
case,  there  is  bound  to  be  some  reflection.  This  reflection  will 
take  place  at  the  level  at  which  the  boundary  condition  is 
applied.  And  so  the  resulting  solution  for  a  non-linear  and  non- 
monochromatic  system  will  be  sensitive  to  the  level  at  which  the 
boundary  condition  is  applied.  Also,  since  the  calculated  wave 
velocity  depends  on  previous  values  of  the  variables,  there  may 
be  a  cumulative  carry  over  of  errors  from  one  time  step  to  the 
next. 

The  boundary  condition  suggested  in  this  paper  is  free  from 
the  above  problems.  However,  we  have  assumed  here  that  the 
numerical  analysis  is  in  terms  of  spectral  Hough  functions.  Hie 
application  of  this  boundary  condition  in  other  methods  of 
solution,  like  finite  difference  methods,  is  a  problem  worthy  of 
further  investigation. 
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Appendix:  Inclusion  of  Coriolis  terms 

Eqs.  (8)  and  (9)  when  Coriolis  terms  are  included  become 

ivu+  -  fv+  +  !£-  =  ft2-  |p  (npQ  f*)  (23) 

ivv+  +  fu+  +  |A1  =  -  |a_  |p(np0  |^)  (24) 

where  u+  =  u  **  yp(  nPc  '  v+  =  v  "  iv  '5T^npo  ^ 
and  f  is  the  Coriolis  parameter. 

The  R.H.S.  of  Eq .  ( 2 '  is  one  order  of  magnitude  less  than 
each  of  the  terms  on  the  L.  H.S.  We  can  therefore  approximate  the 
R.H.S.  of  the  equation  to  be  F^(+’)  where  F^f^’)  is  the 

fg2  3  f  a  vi 

expression  for  jp(.np0  ypj  in  the  solution  of  the  equations 

ivu  -  fv  +  yS—  =  0 
ivv  +  f  u  +  =  0 


This  approximation  in  effect  neglects  products  of  the  diffusion 
coefficients  which  is  in  keeping  with  the  assumption  (in  the  main 
text)  of  small  coefficients  of  diffusion. 

Similarly,  the  R.H.S  of  Eq.  (24)  becomes  F2(<i').  On 
combination  we  have 


9u.  *v+  £(♦’) 

7T  +  7 y  *  Tv  *  0 


(25) 


where 


l -hi 


,  r  -  F,(*')|  +  h  <3f-  '*<♦•>) 


1-  f2/v: 


♦  Ty  t 


,  .ijjr  -  F2<»’>)  ~  I,  Ilf-  *  Fl<*'  >  1 


1  »  f2/\ 


From  the  definition  of  u+  and  v+  Eq.  (25)  becomes 


|iv  -  92  fp  (»<•„  Jp)K|f  *  |i)  +  £(♦•)  -  0 


which  is  the  equation  corresponding  to  Eq,  (13)  with 
£(¥>')  replacing  V2\|»'  , 

Hie  replacement  of  72^'  by  F(^')  is  possible  if  £(♦')  does 
not  contain  operations  like  -jp.  .  However,  in  view  of  the  fact 
that  the  j-p  terms  in  the  F  operator  are  small,  the  replacement 
will  introduce  only  a  very  small  error  which  will  not 
significantly  modify  the  conclusions  we  arrived  at  in  the  main 
text. 
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